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The  irruption  of  gas  and  oil  Into  the  Gulf  of  Mexico  during  the 
Deepwater  Horizon  event  fed  a  deep  sea  bacterial  bloom  that  con¬ 
sumed  hydrocarbons  In  the  affected  waters,  formed  a  regional  ox¬ 
ygen  anomaly,  and  altered  the  microbiology  of  the  region.  In  this 
work,  we  develop  a  coupled  physical-metabolic  model  to  assess  the 
impact  of  mixing  processes  on  these  deep  ocean  bacterial  commu¬ 
nities  and  their  capacity  for  hydrocarbon  and  oxygen  use.  We  find 
that  observed  biodegradation  patterns  are  well-described  by  expo¬ 
nent!  a  i  growth  of  bacteria  from  seed  populations  present  at  low 
abundance  and  that  current  oscillation  and  mixing  processes  played 
a  critical  role  in  distributing  hydrocarbons  and  associated  bacterial 
blooms  within  the  northeast  Gulf  of  Mexico.  Mixing  processes  also 
accelerated  hydrocarbon  degradation  through  an  autoinoculation 
effect  where  water  masses.  In  which  the  hydrocarbon  Irruption 
had  caused  blooms,  later  returned  to  the  spill  site  with  hydrocar¬ 
bon-degrading  bacteria  persisting  at  elevated  abundance.  Interest¬ 
ingly,  although  the  Initial  irruption  of  hydrocarbons  fed  successive 
blooms  of  different  bacterial  types,  subsequent  irruptions  promoted 
consistency  In  the  structure  of  the  bacterial  community.  These 
results  highlight  an  Impact  of  mixing  and  circulation  processes  on 
biodegradation  activity  of  bacteria  during  the  Deepwater  Horizon 
event  and  suggest  an  important  role  for  mixing  processes  in  the 
microbial  ecology  of  deep  ocean  environments. 

oil  spill  |  well  blowout  |  Intrusion  layers 

O  il  and  gas  from  the  Macondo  well  flowed  freely  into  the  deep 
Gulf  of  Mexico  for  a  period  of  83  d  after  the  explosion  and 
sinking  of  the  Deepwater  Horizon  (DWH)  mobile  offshore  drilling 
unit.  The  environmental  release  of  crude  oil  occurred  at  ~\,5  km 
water  depth  with  an  estimated  magnitude  of  4.1  x  106  barrels  (1). 
Large  volumes  of  gas  also  emanated  from  the  ruptured  well,  with 
reported  ratios  of  gas  to  oil  ranging  from  1,600  to  2,800  standard 
cubic  feet  of  gas  (15.6  °C,  I  bar)  per  barrel  of  oil  (2-4).  Mass  fluxes 
estimated  from  these  values  arc  5.4  x  101 1  g  for  liquid  oil  and  1.8- 
3.1  x  101 1  g  for  natural  gases,  defined  here  as  alkanes  with  one  to 
five  carbons. 

Oil  and  gas  entered  the  ocean  initially  through  multiple  open¬ 
ings  in  the  ruptured  riser  pipe  and  later,  from  the  top  of  the 
blowout  preventer  after  the  riser  pipe  was  cut  away  on  June  1, 2010 
(1).  The  hydrocarbon  droplets  ejected  ranged  in  size  from  several 
millimeters  down  to  small  droplets  with  slow  ascent  rates;  771,000 
gal  dispersant  were  applied  at  depth  to  promote  formation  of  such 
small,  slow-rising  droplets.  On  release,  the  bulk  of  oil  and  gas 
began  a  rapid  ascent  from  the  sea  floor,  entraining  sea  water  as  it 
rose.  The  entrainment  of  sea  water  cooled  the  oil  and  gas  rapidly 
(3)  and  initiated  both  dissolution  of  the  soluble  components  and 
formation  of  gas  hydrate.  Kinetically  controlled  chemical  frac¬ 
tionation  seems  to  have  persisted  for  several  hundred  meters  of 
ascent  until  the  entrained  waters  separated  from  the  ascending  oil 
(3,  5,  6).  These  waters  formed  intrusion  layers  between  800  and 
1300  m  water  depth  that  included  dissolved  hydrocarbons,  small 
oil  droplets,  and  dispersant  applied  subsea;  metastable  gas  hy¬ 
drate,  carbon  dioxide,  precipitated  waxes,  and  polar  petroleum 


compounds  arc  also  likely  to  have  been  present  (2,  3,  7-11).  A 
variety  of  processes  (c.g^  slow  ascent  of  small  oil  droplets,  descent 
of  wax  or  other  high  molecular  weight  petroleum  components,  and 
dissociation  of  metastable  gas  hydrate)  likely  altered  plume  com¬ 
position  near  the  source,  leaving  intrusions  dominated  by  the  most 
soluble  compounds,  such  as  gases  (2^4,  9,  10,  12)  and  mono 
aromatic  compounds  (3,  7,  8,  13),  in  addition  to  most  of  the  dis¬ 
persant  added  subsea  (11,  14)  and  a  poorly  defined  quantity  of 
entrained  oil  droplets  (15).  These  intrusion  layer  hydrocarbons 
were  substrates  for  bacterial  respiration,  feeding  a  bacterial  bloom 
that  persisted  for  several  months  (2, 10,  13). 

The  goal  of  this  work  is  to  determine  how  currents  and  mixing 
processes  impacted  the  course  of  hydrocarbon  biodegradation  and 
microbial  ecology  in  the  deep  Gulf  of  Mexico  during  and  after  the 
Macondo  well  blowout.  We  use  a  regional  physical  oceanographic 
model  that  targets  the  motion  of  discrete  water  parcels  in  the  deep 
intrusion  layer.  The  model  was  forced  with  time  course  input  for  26 
abundant  hydrocarbon  compounds  or  classes  and  seeded  with  52 
metabolically  defined  bacterial  types  representing  natural  pop¬ 
ulations  of  both  primary  hydrocarbon  degraders  and  a  secondary 
population  consuming  excreted  intermediary  compounds.  Results 
were  analyzed  with  a  mixing  diagnostic  based  on  the  notion  of 
mesocllipticity  and  mesohypcrbolicity  on  finite  flow  evolution 
timescales  (16).  The  spatial  and  temporal  patterns  of  biode¬ 
gradation,  microbial  community  structure,  and  oxygen  decline  that 
we  obtain  clearly  show  modulation  by  physical  processes  readily 
diagnosed  by  mesocllipticity  and  mesohyperbolicity  and  reveal 
important  effects,  such  as  autoinoculation,  caused  by  recirculation 
of  contaminated  waters  and  their  microbial  populations  about  the 
ruptured  well.  These  results  may  reconcile  disparate  observations 
of  the  physical  dynamics  and  microbial  community  structure  of  the 
deep  plume. 

Model  Development 

We  developed  a  2D  (plus  time)  model  to  combine  horizontal 
circulation  and  mixing  in  the  deep  Gulf  of  Mexico  waters  with 
hydrocarbon  input  that  feeds  bacterial  metabolism  and  cellular 
growth.  We  focused  entirely  on  the  deep  plume  horizon  spanning 
1,000-1,300  m  water  depth,  applying  diagnostics  for  mixing  (16) 
that  have  not  been  used  in  previous  models  of  deep  water  circu¬ 
lation  (17).  We  defined  a  daily  input  flux  of  hydrocarbons  dis¬ 
tributed  evenly  across  depths  in  this  horizon,  with  a  Gaussian 
distribution  about  the  wellhead;  bacterial  populations  of  defined 
initial  abundance,  growth  rate,  and  substrate  preference  were 


Author  contributions;  D.LV,  end  l.M  designed  research;  D.LV„  LM..  S.M..  N  C-Z,  S.L 
and  P.J.H.  performed  research;  D.LV.,  I.M„  VAf,  and  S.L  contributed  new  analytic 
tooh;  DX.V..  LM„  S.M.,  N.C-Z,  S.L  and  PJ.H.  analyzed  data;  and  O.LV.  and  IM  wrote 
the  paper 

The  authors  declare  no  conflict  of  interest 
This  article  is  a  PNAS  Direct  Submission. 

*To  whom  correspondence  should  be  addressed  E-mail:  velentine^geoi.ucsb.adu 
This  article  contains  supporting  information  online  at  www.pna*  orglookup/suppl/doHO. 

107T/pnas  1 108820 lOfVVDC  Supplement  a  I 


20286-20291  |  PNAS  |  December  1 1, 2012  |  vol.  109  |  no.  50 


www.pnas.orgfcgi/doifi  0. 1 073/pnas.  1 108820109 


allowed  to  grow  from  these  inputs.  The  model  tracked  horizontal 
circulation  and  mixing,  bacterial  abundance  and  metabolic  rate, 
and  changes  in  oxygen  and  hydrocarbon  concentrations  resulting 
from  metabolism.  A  general  description  of  the  model  inputs  fol¬ 
lows,  with  detailed  descriptions  provided  in  SJ  Text . 

To  estimate  the  input  of  hydrocarbons  from  the  ruptured 
Macondo  well  into  the  deep  plume  horizon  of  1,000-1300  m  water 


depth,  several  available  sources  were  used.  We  began  with  daily 
flux  estimates  of  oil  and  natural  gas  from  the  well  based  on  ref.  1 
and  subtracted  the  daily  recovery  volumes  reported  from  active 
collection.  We  then  identified  22  compounds  that  were  observed 
or  inferred  (from  solubility)  to  have  partitioned  into  the  deep 
plumes,  and  we  estimated  each  compound’s  daily  emission  from  its 
reservoir  abundance  (3,  4).  For  volatile  compounds,  the  possible 


Fig.  1.  An  analysis  of  the  autolnoculatlon  effect  through  comparison  of  two  water  parcels  with  single-  (red)  and  double-  (blue)  exposure  histories.  (A) 
Trajectories  of  two  water  parcels  for  1  SO  d  starting  April  23,  2010.  (6)  Time  course  of  cumulative  hydrocarbon  Input  to  the  parcels  at  the  proportions  shown  in 
Table  SI .  (O  Time  course  of  hydrocarbon  flux  Into  the  two  parcels  at  the  proportions  shown  In  Table  SI  (D)  Time  course  of  dissolved  oxygen  concentration  In 
the  parcels  attributed  to  hydrocarbon  respiration.  (£)  Time  course  of  respiration  rate  linked  to  methane  consumption  In  the  two  parcels.  (F)  Time  course  of 
bacterial  growth  for  organisms  consuming  methane  (Met  and  Met').  (G)  Time  course  of  respiration  rate  linked  to  consumption  of  25  nonmethane  hydro¬ 
carbons  In  the  two  parcels.  (W)  Time  course  of  bacterial  growth  for  organisms  consuming  25  nonmethane  hydrocarbons  (exdudes  Met  and  MetT  Units  for 
bacteria!  abundance  are  shorthand  for  micromoles  carbon  per  liter. 
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range  of  gas  content  led  to  both  low  (3)  and  high  (2)  assumptions; 
results  presented  derive  from  the  latter.  We  then  applied  partition 
coefficients  as  described  in  SI  Text  to  estimate  the  fraction  of  each 
compound  retained  in  the  deep  plume  and  thereby,  calculated 
a  compound-specific  daily  input  flux  to  the  deep  plume  for  the  full 
duration  of  the  spill.  We  used  a  similar  approach  to  estimate  the 
daily  input  flux  for  the  less-soluble,  longer-chain  hydrocarbons, 
which  were  divided  into  four  classes  based  on  chain  length  and 
assumed  to  have  been  distributed  into  the  deep  plume  as  small 
droplets.  Together,  the  daily  input  fluxes  of  these  26  compounds 
and  classes  arc  the  only  carbon  inputs  to  the  model. 

We  assume  that  each  compound  or  class  is  consumed  by  a  single 
bacterial  lineage,  which  we  define  as  an  operational  metabolic  type 
(OMT),  yielding  26  OMTs  of  primary  consumers  (Table  SI).  We 
assume  the  metabolism  of  primary  coasumers  partitions  carbon 
equally  among  CO2,  biomass,  and  intermediate  metabolites  (in¬ 
cluding  surfactants)  that  escape  the  cell  and  feed  a  secondary 
microbial  community.  The  escape  of  intermediate  metabolites  is 
typical  of  hydrocarbon  degrade rs  (18-20),  presumably  because 
such  intermediates  are  hydrophobic  and  readily  cross  the  cell’s 
membrane.  Secondary  consumers  arc  uniquely  linked  to  the  in¬ 
termediate  produced  by  the  associated  primary  consumer,  yielding 


26  OMTs  of  secondary  consumers  (Table  SI).  These  secondary 
consumers  are  assumed  to  partition  their  substrate  evenly  into 
CO2  and  cellular  biomass.  Biomass  generated  by  both  primary  and 
secondary  consumers  is  assumed  to  be  remincralizcd  according  to 
a  first-order  decay  process  (k  =  2  x  1  (T3  IT1)  when  substrate 
concentration  no  longer  supports  growth.  Each  OMT  is  assumed 
to  be  present  throughout  Gulf  of  Mexico  deep  waters  at  an 
abundance  of  104  viable  cells  LT1  (0.1%  of  the  total  microbial 
community  of  107  cells  L“*)  (13),  except  for  methanotrophs,  which 
were  assumed  to  be  present  at  103  viable  cells  L“,.The  background 
community  was  assumed  to  be  present  at  0.1  pmd-C  (107  cells 
L-1  at  120  fg  C  per  cell).  All  OMTs  mapped  to  one  of  the  following 
bacterial  famflies^genera  (substrate  preference  in  brackets):  Methyl - 
ococcQccac  [methane  (2,  21)],  Oceanospirillales  [alkanes  (13)], 
Cycloclasticus  [aromatics  (18)],  McthylophagalMethyiophilaceac 
[methane-derived  intermediates  (2)],  and  Cotwcllia  [propane 
ethane,  and  higher  alkane-derived  intermediates  (21, 22)].  Where 
possible,  OMTs  were  mapped  to  bacterial  species,  and  published 
data  on  these  species  were  used  to  estimate  doubling  times  and 
cell  masses.  Each  OMT  was  assumed  to  grow  with  a  characteristic 
doubling  time  after  it  encountered  the  appropriate  substrate, 
limited  only  by  substrate  depiction.  The  estimated  C  masses  and 


Fig.  2.  Comparison  of  mkrobial  community  dynamics  in  water  parcels  with  different  exposure  histories,  all  starting  April  23,  2010.  (A)  Time  course  change  in 
hydrocarbon  flux  (red)  and  respiration  rate  (blue)  shown  in  Top,  dissolved  oxygen  concentration  (red)  and  bacterial  abundance  (blue)  shown  In  Middle,  and 
relative  composition  of  the  microbial  community  shown  In  Bottom  for  a  parcel  experiencing  a  single  exposure.  (£)  Time  course  changes  as  described  for  A  for 
a  parcel  experiencing  a  triple  exposure.  (O  Time  course  changes  as  described  in  A  for  a  parcel  experiencing  a  double  exposure.  (D)  Tabulated  legend 
identifying  the  OMTs,  their  putative  phylogenetic  affiliation,  and  their  ecosystem  function. 
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Fig.  3.  impact  of  recirculation  on  hydrocarbon  abundance,  bacterial  popuiatioa  and  metabolism  (A  and  B )  Spatial  distribution  of  bacterial  abundance  before 
CA)  and  during  ( B)  the  early  stages  of  an  autoinocuiation  event.  (O  Time  course  change  In  average  abundance  for  bacteria  consuming  nonmethane  hydro¬ 
carbons  and  the  average  summed  concentration  of  these  chemicals  integrated  over  the  computational  domain.  (0  and  E)  Spatial  distribution  of  hydrocarbon 
respiration  rate  before  (0)  and  during  (E )  the  early  stages  of  the  same  autoinocuiation  event  (f)  Time  course  change  in  the  average  respiration  rate  for  bacteria 
consuming  nonmethane  hydrocarbons,  integrated  over  the  computational  domain.  Because  the  size  of  the  computation  domain  Is  2#  latitude  x  2*  longitude,  a 
small  number  of  parcels  exit  the  domain  near  the  end  of  the  simulation. 


doubling  times  for  each  of  the  52  OMTs  are  provided  in  SI  Text 
along  with  the  associated  reasoning. 

Results  and  Discussion 

The  model  predicts  that  the  2  x  101 1  g  hydrocarbons  injected  into 
the  deep  Gulf  of  Mexico  drove  total  bacterial  productivity  of  101 1  g 
C  (-Kr3  cells).  Such  productivity  is  sufficient  to  double  the  bac¬ 
terial  population  in  the  deep  plume  horizon  ( 1,000-1300  m  depth) 
over  an  area  of  30,000  km2  and  generate  a  regional  oxygen 
anomaly  of  I012  g.  Published  observations  suggest  that  this  process 
began  within  weeks  of  the  initial  irruption  (9,  10,  13)  and  pro¬ 
gressed  over  a  period  of  -120  d,  by  which  point  the  natural  gases 
had  been  consumed  and  the  fluorescence  anomaly  in  the  deep 
plume  had  been  strongly  attenuated  (2).  These  observations,  made 
during  the  DWH  event,  serve  to  ground  truth  our  models’  pre¬ 
dictions  of  the  temporal  and  spatial  distributions  for  microbial 
growth  and  metabolism. 

Under  our  model,  when  a  water  parcel  first  encounters  a  hy¬ 
drocarbon  source,  a  seed  population  of  active  hydrocarbon- 
degrading  cells,  present  at  low  abundance  in  all  waters,  begins  to 
grow.  As  growth  progresses,  OMTs  outpace  the  influx  of  their 
preferred  substrates,  giving  sequential  pulses  of  respiration  as 
different  OMTs  flourish  and  then  fade  (Fig.  1).  The  timing  of 
respiratory  pulses  is  controlled  primarily  by  the  growth  rates  of  the 
responsible  OMTs,  initial  abundance  of  active  cells,  and  duration 
of  source  exposure.  The  model  shows  initial  pulses  of  growth  and 
respiration  linked  to  butane  and  pentane  followed  by  sequential 
pulses  from  consumption  of  propane  and  other  short-chain  alka¬ 
nes,  ethane,  aromatic  hydrocarbons,  and  finally,  methane. 

By  injecting  the  hydrocarbons  into  moving  parcels  of  water,  the 
model  generates  a  complex  spatial  pattern  of  bacterial  growth  and 
metabolism  (Movie  SI)  and  reveals  associated  feedback  mecha¬ 
nisms.  Tracking  individual  water  parcels  as  they  advcct  reveals  that 
many  parcels  come  into  contact  with  the  hydrocarbon  source  on 
multiple  occasions  (Fig.  1).  A  parcel’s  first  encounter  with  hydro¬ 
carbons  feeds  staggered  pulses  of  respiration  and  successive 
blooms  of  individual  OMTs  (Fig.  2A)y  and  residual  hydrocarbon¬ 


degrading  OMT  abundances  remain  elevated  above  background. 
We  apply  the  term  autoinocuiation  to  describe  the  priming  effect 
that  occurs  when  hydrocarbons  are  introduced  into  such  a  pre¬ 
viously  exposed  parcel.  When  hydrocarbons  irrupt  into  the  pre¬ 
viously  exposed  parcel,  the  elevated  abundance  of  hydrocarbon- 
degrading  OMTs  allows  simultaneous  consumption  of  different 
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Fig.  4.  Comparison  of  physical  and  biological  features.  (A)  Spatial  distri¬ 
bution  of  bacterial  abundance  on  June  13,  2010.  ( B)  Hypergraph  map  on 
June  13,  2010.  (O  Spatiai  distribution  of  bacterial  abundance  on  June  30, 
2010.  (0)  Hypergraph  map  on  June  30,  2010  For  the  hypergraphs,  blue 
indicates  meso hyperbolic  regions,  red  Indicates  mesohyperboikity  with 
shear,  and  green  indicates  elliptic  regions  in  which  eddies  are  formed.  The 
scaied  units  for  the  hypergraphs  are  defined  in  ref.  16. 
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compounds  at  rates  substantially  faster  than  occurred  after  initial 
exposure.  This  effect  is  apparent  in  Fig.  1  by  comparing  bacterial 
abundance  (Fig.  1 H)  with  respiration  pulses  (Fig.  1G)  for  single 
(red)  vs.  double  (blue)  exposure  cases,  noting  that,  after  the  second 
exposure  occurs  at  53  d,  previously  irregular  growth  and  re¬ 
spiratory  pulses  are  collapsed  into  single  peaks.  In  the  case  of 
methane,  growth  (Fig.  IE)  and  respiration  (Fig.  IF)  associated 
with  the  second  exposure  are  subsumed  into  a  single  pulse  paced 
by  the  timing  of  the  initial  exposure.  Collapse  of  the  respiratory 
succession  leads  to  subtle  changes  in  bacterial  community  struc¬ 
ture  (Fig.  2),  thereby  avoiding  the  dramatic  changes  in  OMT 
abundance  and  community  structure  that  followed  the  initial  ex¬ 
posure.  Autoinoculation  seems  to  have  played  an  important  role  in 
suppressing  hydrocarbon  buildup  in  the  deep  plume;  because 
discharge  from  the  well  continued  for  so  long,  roughly  one-half  of 
the  total  irruptive  hydrocarbon  discharge  entered  preexposed 
waters.  In  one  instance  (Fig.  3),  the  model  predicts  that  waters 
returning  from  the  northwest  in  mid-June  persisted  over  the 
wellhead  until  mid-July  when  discharge  ceased,  resulting  in  what 
we  term  an  autoinoculation  event.  Crucially,  the  accelerated  res¬ 
piration  resulting  from  this  event  both  drove  down  levels  of  non¬ 
methane  hydrocarbons  and  stabilized  the  overall  respiration  rate, 
which  had  previously  been  more  variable  (Fig.  3). 

Currents  and  mixing  processes  modulated  the  distribution  and 
activity  of  bacteria  in  the  deep  Gulf  of  Mexico,  frequently  rein¬ 
troducing  water  parcels  to  the  source  area.  Analysis  of  our  model 
simulations  shows  that  mixing  features  of  the  flow  field  influenced 
the  chemical  and  microbial  distributions  at  plume  depth.  Com¬ 
parison  of  figures  representing  bacterial  abundance  and  hypergraph 
maps  (Fig.  4  A  and  C  vs.  Fig.  4  B  and  D)  indicates  that  blue  mes- 
ohypcrbolic  features  in  the  hypergraph  map  accumulate  high  bac¬ 
terial  density;  such  features  act  as  convergence  zones  for  the 
advccted  reacting  field.  The  green  mcsoelliptic  zones  in  the  flow 
field  cause  a  swirling,  recurrent  motion  of  the  water  and  the 


associated  bacterial  population  (Fig.  4  A  and  £).  This  association 
indicates  that  the  spatial  distribution  of  the  nonlinearly  reacting 
flow  field  is  strongly  affected  by  the  velocity  field  mixing  features. 
Their  importance  for  the  autoinoculation  mechanism  is  particularly 
evident  in  Fig.  4/J,  where  the  hypergraph  map  indicates  large  mixing 
events  around  the  well  that  arc  visualized  by  the  mixture  of  red  and 
blue.  Mcsoelliptic  zones  influenced  bacterial  communities,  such  as 
in  late  May  of  2010,  when  such  a  zone  promoted  a  buildup  of 
hydrocarbons  in  the  vicinity  of  the  wellhead  (Fig.  5).  This  recircu¬ 
lation  feature  may  have  provided  a  false  impression  of  a  plume 
emanating  to  the  southwest,  which  was  reported  for  that  time  (13). 
Mcsohyperbolic  features  may  have  also  been  inadvertently  identi¬ 
fied  during  the  spill,  such  as  in  June  of  2010  (Fig.  SI),  when  a  sin¬ 
uous  plume  was  tracked  to  the  west  of  the  wellhead  (7). 

Wc  find  good  agreement  between  the  model  and  available  data 
with  respect  to  both  plume  motion  and  location  and  microbial 
community  distribution  and  activity.  The  plume  was  reported  to 
emanate  to  the  southwest  (2, 7-9)  at  some  times  and  the  northeast 
at  other  times  (10),  eventually  migrating  to  the  southwest  (2)  in 
August  of  2010;  the  modeled  current  oscillates  between  the 
northeast  and  southwest,  with  eventual  migration  to  the  southwest 
in  August  of  2010.  Samples  collected  close  to  the  ruptured  well  in 
June  of  2010  showed  different  degrees  of  hydrocarbon  degradation 
(7,  9,  10,  13);  the  model  predicts  that  previously  exposed  water 
returned  to  the  wellhead  at  this  time  (Fig.  3),  leading  to  mixed 
levels  of  biodegradation  in  the  vicinity  of  the  wellhead.  Samples  of 
these  waters  revealed  blooms  of  Cotwctiia  and  a  biodegradation 
preference  for  propane  over  ethane  over  methane  (9,  10);  the 
model  predicts  concurrent  blooms  of  OMTfc  mapped  to  Cohicltia , 
which  have  been  linked  to  ooasumption  of  ethane  and  propane 
(21).  The  model  provides  for  abundant  methanotroph  and  meth- 
yiotroph  production,  which  was  observed  in  September  and  Octo¬ 
ber  of  2010  (2),  and  it  independently  arrives  at  the  debated  scenario 
(23-25)  in  which  methanotrophs  and  methylotrophs  constitute 


Fig.  5.  Spatial  and  temporal  context  for  development  of  bacterial  communities  reported  by  Hazen  et  ai.  (13)  for  May  30,  2010.  (A)  Spatial  distribution  of 
bacterial  abundance  for  primary  hydrocarbon  consumers  on  May  30, 2010.  ( B )  Spatial  distribution  of  the  dissolved  oxygen  anomaly  formed  from  hydrocarbon 
respiration  on  May  30.  2010.  (C)  Spatial  distribution  of  totai  hydrocarbons  (sum  of  26  compounds  or  dasses)  on  May  30,  2010.  (D)  Spatial  distribution  of 
respiration  rate  for  hydrocarbon  oxidation  on  May  30,  2010.  (£)  Seven-day  hypergraph  for  May  30,  2010,  showing  mixing  regimens  and  highlighting 
a  mesoeliiptlc  zone  where  an  eddy  has  set  up  to  the  southwest  of  the  wellhead.  Scaled  units  are  defined  In  ref.  16.  (F)  One  hundred  fifty-day  trajectory  of 
a  water  parcel  located  near  the  position  of  sampling  (1 3)  on  May  30,  2010.  (G)  Time  course  changes  for  the  water  parcel  for  which  the  trajectory  is  shown  In  F 
(as  per  Fig.  2);  May  30,  2010  Is  shown  as  a  dashed  line.  The  location  of  the  Macondo  well  Is  denoted  by  a  black  drde  In  A-£.  Samples  collected  by  Hazen  et  al. 
(13)  on  this  day  were  within  the  symbol  denoting  the  wellhead.  Note  the  iog  scales  for  A,  B.  and  D. 
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a  significant  portion  of  the  bacterial  community,  even  with  no  de¬ 
tectable  oxygen  anomaly  (c.g.,  respiration  of  0.05  pM  CR4  yields 
a  maximum  oxygen  anomaly  of  0.1  pM,  below  the  detection  limit; 
the  resulting  bacterial  biomass,  a  maximum  of  0.025  pmol  C  L”1, 
would  constitute  up  to  20%  of  the  background  community  with  no 
apparent  oxygen  anomaly,  coasistcnt  with  our  previous  observa¬ 
tion)  (2).  Hazen  et  al.  (13)  observed  that,  on  May  30,  90-95%  of 
the  microbial  community  in  the  peak  plume  horizon  southwest 
of  the  wellhead  was  OceanospiriUalcs\  the  model  predicts  ~80%  of 
the  bacterial  population  as  OMTs  mapped  to  Occanospirillalcs  at 
this  place  and  time  (Fig.  5).  Likewise,  column-integrated  bacterial 
productivity  in  the  deep  plume  at  this  place  and  time  was  modeled 
to  be  0.16  mol-C  nT2  (ranging  from  0.05  to  050  mol-C  m~2  within 
5  km);  this  finding  is  in  good  agreement  with  the  calculated  value  of 
0.20  mol-C  m“2 * *  based  on  a  3.9  x  107 8 * 10 * 12*cclls  L'1  deasity  of  Ocean - 
ospirillales  (13),  cell  carbon  content  (Table  SI)  of  5.0  x  10" 13  g-C 
cell”1,  and  local  plume  thickness  of  120  m  (13).  The  locations  of 
observed  oxygen  anomalies  were  also  consistent  with  model  results 
as  in  Movie  S2.  In  this  simulation,  we  compare  the  presence  or 
absence  of  oxygen  anomalies  from  190  hydrocasts  conducted  during 
the  time  of  active  discharge  (2,  7-10,  13,  26)  with  the  contempo¬ 
raneous  model  output.  We  find  a  72%  match  rate;  kilometer-sealed 
heterogeneity  may  account  for  up  to  93%  of  the  remaining  obser¬ 
vations  for  which  the  model  incorrectly  predicted  as  anomalies. 

Nonetheless,  some  discrepancies  arc  apparent  between  model 
and  observation,  particularly  as  the  time  course  modeled  grows 
longer.  Most  notably,  the  model  undcrprcdicts  the  velocity  of 
plume  migration  in  August  and  September  of  2010,  with  a  con¬ 
current  latitudinal  offset;  although  the  model  properly  moves  the 
plume  to  the  southwest,  observations  show  the  plume  arriving  ~2 
wk  earlier  and  ~20  km  farther  south  than  the  model’s  prediction. 
The  model  also  fails  to  capture  the  small-scale  spatial  heteroge¬ 
neity  apparent  from  response  efforts.  This  discrepancy  is  attrib¬ 
utable  to  the  model’s  ~4-km  grid  size,  the  Gaussian  input  of 
hydrocarbons,  and  the  even  distribution  of  hydrocarbons  within 
the  depth  horizon.  Unfortunately,  because  sampling  typically  tar¬ 
geted  the  anomalies,  local  heterogeneity  also  complicates  direct 
comparison  of  model  output  and  chemical  data.  Although  model 
tuning  improves  our  ability  to  match  observations,  we  find  gTeat 
use  in  explaining  biogcochcmical  and  microbiological  observations 


from  fundamental  principles  of  chemistry,  microbiology,  and 
ocean  physics  without  the  need  for  tuning  or  optimization. 

In  light  of  the  model’s  overall  success,  we  considered  how  it 
might  inform  the  interpretation  of  results  from  previous  studies. 
First,  the  inclusion  of  physical  dynamics  corroborates  a  previous 
hypothesis  (23,  25)  for  why  the  extensive  and  persistent  oxygen 
depletion  predicted  in  the  deep  plume  horizon  by  Joye  et  al.  (9) 
was  not  observed.  Second,  the  physical  dynamics  also  suggest  a 
more  complex  flow  field  than  the  laminar  flow  used  by  Hazen 
et  al.  (13)  to  estimate  hydrocarbon  half-lives  associated  with  the 
bloom,  thus  calling  the  accuracy  of  the  reported  half-lives  into 
question.  Finally,  the  modeled  flow  regimen  predicts  that  a  sinu¬ 
ous  plume  similar  to  the  plume  identified  in  the  work  by  Camilli 
et  al.  (7)  may  not  have  come  directly  from  the  wellhead  but  rather, 
had  first  traveled  northeast  before  recirculation  to  the  southwest 
(Fig.  SI).  Although  its  predictions  cannot  be  conclusive,  this 
model,  thus,  provides  a  useful  holistic  context  for  the  inter¬ 
pretation  of  biological,  chemical,  and  physical  observations. 

Our  results  suggest  that  currents  and  mixing  processes  acceler¬ 
ated  hydrocarbon  biodegradation  and  structured  microbial  ecology 
during  the  DWH  event  through  priming  of  the  native  microbiota 
residing  in  water  parcels  that  recirculated  around  the  ruptured  well. 
Although  memory  effects  are  common  in  enclosed  environments 
and  physical  processes  are  known  to  shape  microbial  biogeography 
(27),  we  suggest  that  such  a  dynamic  autoinoculation  effect  may 
be  of  broad  importance  to  deep  ocean  environments.  The  impacts 
on  ocean  biogcochemistiy  and  microbiology  may  encompass  point 
sources,  such  as  hydrothermal  vents  and  hydrocarbon  seeps,  but 
may  also  include  suboxic  water  masses,  organic-rich  sediment 
deposits,  flushing  of  silled  basins,  and  regions  receiving  high-export 
flux.  For  ocean  waters  susceptible  to  contamination  from  events 
such  as  the  DWH,  a  major  challenge  remains  to  predict  how 
physical  dynamics,  chemistry,  and  microbiology  interact  to  set  the 
timescale  and  extent  of  the  microbial  response. 
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